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基于 CEI 测 量 的 GEO 目 标 快速 轨道 恢复 
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\ 算 法 
法 . 利用 亚太 
30 min 短 弧 定 轨 和 10 min 预 报 
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连 线 干涉 测量 (Connected Element Interferometry,，CEI) 是 一 种 全 
标的 跟踪 监视 . 地 球 静 止 轨道 (Geostationary Earth Orbit, GEO) E.A 


L 动 后 的 快速 轨道 恢复 能 力 对 于 监视 预警 极为 
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Ec C, 在 附加 


号 卫星 的 CEI 仿 


先 验 轨道 约束 的 短 弧 定 轨 直 
真 和 实测 数据 进行 了 短 弧 定 轨 和 


的 卫星 
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1 引言 

地 球 静 止 轨道 (Geostationary Earth Orbit, 
GEO) 因 其 轨道 资源 的 独特 性 与 稀缺 性 , 成 为 各 国 
相互 争夺 的 重要 太空 资源 上 4. 由 于 摄 动 力 等 因素 
的 影响 GEO 目 标 需 要 频繁 地 实施 机 动 以 保持 在 
分 配 的 定点 区 间 内 . 这 种 位 保 机 动 破坏 了 常规 定 轨 


精度 的 连续 一 致 性 ， 


因此 机 动 后 的 快 


GEO 目 标的 常规 管理 至 关 重 要 回 ， 


针对 非 合作 
随机 的 特殊 任务 机 动 ， 因 


GEOH 


标的 监视 预警 , 还 要 应 对 各 
此 监测 系统 的 快速 响 


速 轨道 恢复 对 


应 能 力也 极为 迫切 . 连 线 干涉 测量 (Connected El- 


ement Interferometry, CEI) 
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GEO 目 标 且 


日 常 机 动 监 视 和 轨道 确定 四. 


短 弧 定 轨 . I 


种 基于 
电信 号 的 被 动 测量 方式 , 可 以 全 天 候 跟 踪 非 合作 
\ 有 隐 珊 性 , 已 被 应 用 于 GEO 卫 星 的 


卫星 下 行 


但 是 ，GEO 目 标 机 动 后 的 快速 轨道 恢复 属于 
于 GEO 目 标的 静 地 特性 
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分 量 精度 均 优 了 
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F4 km, 


报 , 实验 结果 表明 , 采 
能 够 满足 非 合 作 GEO 目 


全 天 候 的 被 动 测 
要 频繁 机 动 以 保 
EH. 针对 基于 CEI 的 GEO 短 弧 定 负 
础 上 ,提出 了 轨道 半 长 有 


NFE KPI AME 刘 隆 迪 


技术 , CH 
持 轨 位 或 完 
L 和 预报 ,分 


初 值 的 自 适应 优化 方 


标 机 动 后 1 


优化 后 的 半 长 轴 初 值 ， 


天 速 轨道 恢 


轨道 力学 约束 弱 , 短 弧 定 轨 通 常 有 较 强 奇异 性 ， 从 


Tft xt EA iul 


12 711380) 
献 [8] 提 


用 C 频 


短 弧 定 轨 和 5 min 预 报 


精度 的 影 
大 于 3 km; 双 站 8 h 弧 长 定 轨 ， 
居 累 积 时 间 长 , 对 测 站 构 型 要 求 较 
于 9 参数 星 历 拟 合 的 短 弧 运动 学 定 轨 方 
法 , 针对 参数 拟 合 的 奇异 问题 提出 了 解决 措施 . 利 
0D 星 的 实验 表明 , 10 min 
的 位 置 精度 均 优 于 20 m, 但 


FE 精度 . 文献 四 依据 测 角 定 轨 原 理 ， 分析 
了 测 轨 弧 长 和 测 站 数 对 非 合 作 GEO 卫 8 


E 定 轨 预 报 


LT 


向 . 小 于 4 h 弧 段 的 单 站 定 轨 , 位 置 误差 均 
精度 可 优 于 3 km. 但 


IH 
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段 测 距 对 北斗 GEOJ 


该 方法 仅 适 


外 
性 , 提 


束 的 短 弧 定 轨 
表现 为 形 亏 , 而 半 长 轴 与 平 经 度 有 较 强 
加 先 验 方差 约束 
并 针对 轨道 : 
定 轨 结果 影响 大 的 现象 , 提出 了 半 长 轴 初 值 的 自 适 


敏感 ， 


相关 性 , 表现 为 数 亏 . 


于 合作 目标 . 
| 对 GEO 非 合作 
出 一 


' 自 适应 优化 半 长 区 


标的 CEI 短 弧 定 轨 的 奇异 
和 附加 先 验 轨道 约 


因此 ， 


法 的 病态 性 ， 


通过 附 


“长 而 


法 . CEI 观 测量 对 卫星 径 向 信息 不 


误差 


由 的 初 值 对 
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应 优化 方法 . 利用 郑州 的 CEI 测 轨 系 统 对 亚太 七 号 
卫星 的 实验 , 验证 了 该 算法 对 短 弧 定 轨 和 预报 的 精 
度 提 升 , 为 GEO 机 动 后 的 快速 轨道 恢复 提供 了 改进 
依据 . 


2 ”CEI 短 弧 定 轨 原理 
2.1 WASE 

CEI 系 统 的 测量 原理 如 图 1 所 示 , CEI 观 测量 是 
基线 两 端 天 线 接收 到 的 下 行 信号 经 互相 关 后 得 到 
的 不 足 一 个 波长 的 相位 差 e-120. 图 中 基线 AB 上 的 
量 测 方程 为 : 


$ +AN = pa — pp = |Ra(t1) — r(to)| 
|i (t2) — r(to)| + eAtacact 
ADatm T A Dins TE , (1) 


式 中 ,% 表 示 相 位 差 观 测 值 (以 长 度 为 单位 )， 和 和 NN 
为 信号 的 波长 和 整 周 模糊 度 ; pA 与 pe 表示 卫星 至 
主 站 A 与 副 站 B 的 距离 . r(to) 表 示 卫 星 信号 发 出 时 
刻 如 卫星 位 置 向 量 , RA 和 Re 则 分 别 是 主 副 站 在 接 
收 时 刻 右 和 如 的 位 置 向 量 . Ataoek 表 示 站 间 钟 同步 
偏差 ; Apatm 和 Apins 为 大 气 传播 延迟 和 设备 延迟 的 
站 间 差 ; = 为 观测 噪声 . 

连 线 式 CEI 通 常 为 数 十 米 至 数 十 干 米 的 短 基 
线 , 相对 于 高 轨 静 地 GEO 卫 星 , 相位 差 观 测量 可 看 


作 是 一 种 ire Hd eil i24. 为 了 使 GEO 
卫星 的 天 球面 指向 精准 , 至 少 需 要 两 条 相交 基线 ， 
AET 正 交 基线 . 


要 指出 的 是 , 使 用 CEI 相 位 差 观测 的 前 提 是 
EAEEREN. 由 于 常规 机 动 对 GEO 轨 位 的 改 
变通 常 在 干 米 级 , 且 基 线 越 短 对 这 种 变 轨 越 不 敏感 ， 
因此 短 基线 条 件 下 , 可 以 采用 机 动 前 的 模糊 度 解 算 
值 . 


2.2 IEAA 

GEO 卫 星 的 初 轨 状态 量 , 在 地 心 惯性 系 (Earth 
Center Inertial Coordinates System, ECI) 下 通常 
采用 专用 于 小 偏心 率 和 小 倾角 的 第 二 类 无 奇 点 根 
数 表示 , Moo = (a, e, à, Ao). 它 与 经 典 开 普 勒 轨道 
根 数 的 关系 如 下 [25713]; 


cu 


学 报 1 期 


工 


e = (es, e). = (ecos ®, esind)? , 


i= (ini) = 


ào = R+w+ Mo, 


式 中 , e 为 二 维 偏心 率 向 量 ; i 是 二 维 倾角 向 量 ; w = 
N 十 w; 和 o 称 为 初始 时 刻 的 平 经 度 ， 从 春分 点 计量 
至 该 时 刻 的 卫星 位 置 ; Mo 为 初始 时 刻 的 平 近 点 角 . 
开 普 勒 轨道 根 数 (a,e,i, 2, w, M) 分 别 表 示 轨 道 半 
长 轴 、 仿 心率、 轨道 倾角 、 升 交点 亦 经 、 近 点 角 
距 和 平 近 点 角 . 


ee 


Radio 


(sinicosQ,sinisinf2)", (2) 


y! Sub Station C 


图 1 CEI 系 统 测量 原理 示意 图 


Fig.1 Schematic diagram of CEI system measurement 


principle 


短 弧 定 轨 的 力 模 型 可 仪 考虑 二 体 问题 , 弧 段 内 
任意 时 刻 t; 的 线性 化 观测 误差 方程 可 表示 为 : 


Uj 一 Ajóco = Aó; ; (3) 
其 中 6oo 是 初 轨 状 态 的 改正 数 向 量 ， A9; 是 观测 值 
与 理论 计算 值 的 差 值 ( 即 0-C), v; 称 为 观测 残 差 ， 系 
数 矩 阵 为 
| (99j Or; Oc; 
A-($S) (m2) (a2) . 的 


1x3 3x6 
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采用 最 小 二 乘 平 差 , 即 满足 V7PV = min, 则 
可 求解 初 轨 状 态 量 的 最 佳 改正 数 向 量 


660 = —(ATPA) ! AT PI, (5) 
其 中 , V 是 观测 残 差 组 成 的 列 向 量 , P 是 观测 权 甜 
阵 , l 是 自由 项 , 即 弧 段 内 Ag9; 组 成 的 列 向 量 ; 由 于 


GEO 对 地 的 构 型 基本 不 变 , 权 和 矩 阵 忆 可 取 单 位 阵 


3 短 弧 定 轨 的 改进 算法 

由 于 精度 要 求 不 高 ,， 短 弧 定 轨 也 可 以 直接 在 
地 心地 固 系 (Earth Centered Earth Fixed System, 
ECEF) 下 实现 . 二 体 问题 下 的 ECEF 初 轨 状 态 , 1X. 
需要 将 第 二 类 无 奇 点 根 数 的 历 元 平 经 度 和 No 改换 为 
历 元 时 刻 的 定点 经 度 和 oo, RU 


Aco 二 Ao So ? (6) 


式 中 So 是 to 时 刻 的 格林 尼 治 时 角 . 
若 忽 上 略 小 倾角 和 小 偏心 率 的 影响 , 任意 时 刻 t 
的 定点 经 度 Xe 能 够 反映 GEO 卫 星 的 对 地 变化 , 即 


(7) 


式 中 和 me 分 别 是 卫星 运动 和 地 球 自 转 的 平均 角 速 
率 , 由 于 GEO 的 静 地 特性 , 二 者 仅 有 微小 偏差 . 
此 外 , 依据 GEO 工 作 卫 星 的 轨 位 分 配 机 制 , 定 
点 经 度 和 通常 规定 在 标 称 经 度 附 近 的 土 0.1° 内 变化 ， 
这 也 为 该 参数 的 先 验 约束 带 来 了 便利 29l. 
因此 ，ECEF 下 的 短 弧 定 轨 初 轨 状 态 量 90 
(a, e, 4, Xeo)T， 相 应 的 平 差 求 解 公式 也 需要 略 做 修 
改 . 


3.1 


Ae = No (n — ne) (t — to), 


病态 性 分 析 
基于 CEI 的 GEO 短 弧 定 轨 有 较 强 的 病态 性 , 表 
现 为 相互 关联 的 形 亏 和 数 亏 . 

(1) 因 CEI 测 角 本 质 导 致 的 形 亏 
CEI 的 短 基线 使 得 高 轨 卫 星 下 行 信号 的 球面 波 
传播 可 理想 化 为 平面 波 传播 . 此 时 ，CEI 相 位 差 表 


ERWE: 基于 CEI 测 量 的 GEO 


标 快速 轨道 恢复 1H 


定 了 其 较 弱 的 轨道 力学 约束 , 短 弧 情 况 下 约束 性 更 
加 弱化 9. 在 缺乏 卫星 测 距 观测 量 的 情况 下 , 仅 基 


mi 


于 CEI 的 GEO 短 弧 定 轨 模型 鲁 棒 性 差 , 导致 无 法 解 
SEC. 


Baseline 


2. CEIF MI 


R] 


ZIN i 


Fig.2 Angular measurement of CEI 


(2) 因 轨道 状态 参数 相关 性 引起 的 数 亏 
由 开 普 勒 第 三 定律 mw?a? = u (4 是 地 球 引力 常 
数 ), 定点 经 度 和 A。 与 半 长 轴 a 的 改正 数 之 间 存 在 线性 


漂移 , BI 


ee snAalt i3. - 748) 


AF, AXe 为 定点 经 度 的 变化 量 , Am 为 卫星 运动 平 
均 角 速率 和 地 球 自转 平均 角 速 率 的 差 值 ， Ao 为 半 
长 轴 的 变化 量 . 
前 面 分 析 了 CEI 的 形 亏 使 得 短 弧 定 轨 无 法 改进 
半 长 轴 a, 这 里 半 长 轴 的 初 值 偏差 又 将 直接 引起 定 
点 经 度 的 线性 漂移 . 因此 , 历 元 时 刻 定点 经 度 和 oo 的 
解 算 精 度 与 半 长 轴 初 值 的 精度 共有 密切 关系 , 若 处 
理 不 好 两 者 均 难 以 得 到 改进 . 


改进 算法 
为 了 削弱 短 弧 定 轨 的 病态 性 , 必须 引进 先 验 轨 


3.2 


现 为 卫星 相对 于 基线 的 方向 变化 信息 CEI 干涉 测 
角 示 意图 如 图 2 所 示 . 因此 ， 这 种 “ 测 角 ”类 型 观测 
量 对 卫星 的 径 向 位 置 不 敏感 ， 且 GEO 静 地 特性 决 
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道 约 束 . 除了 增加 常规 的 参数 先 验 协 方差 约束 , 还 
必须 确定 最 优 的 半 长 轴 初 值 . 

(1) 附 加 先 验 约束 
依据 GEO 的 机 动 特性 和 定点 控制 区 规定 ， 
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初 轨 参 数 的 先 验 均 方差 约束 为 : 半 长 轴 a 取 土 20 km, 

历 元 时 刻 的 定点 经 度 和 oo 和 倾角 分 量 i 取 土 0.1°. 
附加 先 验 约束 的 初 轨 改 正法 方程 为 : 


ôo = -(A'PA* P3, ) A'PI, (9) 
式 中 ，6G0 为 附加 先 验 约束 的 ECEF 下 短 弧 定 轨 初 
轨 状 态 量 的 最 佳 改正 数 向 量 , Pj 为 ECEF 下 短 弧 
定 轨 初 轨 状 态 量 的 先 验 权 和 矩阵 ; 其 中 ， 若 取 单位 权 


方差 为 1, 则 先 验 权 和 矩阵 Pj 的 对 角 线 元 素 为 : 


diag (Pa) = 
{(20km)? 0 0 0 (0.1°) 


(oe 
法 矩阵 的 正定 


(10) 


显然 , 先 验 约束 的 实质 就 是 增强 
TE, 从 而 改善 矩阵 求 逆 的 病态 问题 . 
(2) 自 适应 优化 半 长 轴 初 值 
常规 定 轨 对 初 值 精度 要 求 不 高 . 以 GEO 
为 例 ， 半 长 轴 的 初 值 可 取 为 静止 轨道 的 标 称 值 
(42165.76 km), 定点 经 度 的 初 值 Mo 可 取 为 定点 经 
度 的 标 称 值 , 其 他 轨道 参数 均 取 0 即 可 . 
但 是 , 在 实际 解 算 中 发 现 , 若 半 长 轴 初 值 在 标 
称 值 附近 一 定 范围 内 随机 取 值 ， 解 算 精 度 会 出 现 
较 大 差异 . 在 没有 其 他 辅助 测量 手段 的 前 提 下 , 本 
文 经 验 性 地 设计 了 一 种 自 适 应 选取 半 长 轴 初 值 的 
优化 算法 . 时 如 下 : 截取 静止 轨道 标 称 值 的 
+20 km 作为 探测 窗 H, 


pE 


-FL 
~ 


T 
z 
z 


流利 Hj 
以 0.2 km 为 步 长 遍历 窗口 


FIFKA 


武 初 值 ， 分 别 解 算 轨 道 . 考察 解 算得 到 
的 偏心 率 序列 ( 即 e = Ve2 + e2), 将 偏心 率 极 小 值 
对 应 的 半 长 轴 作 为 半 长 轴 的 最 佳 初 值 . 偏心 率 极 小 
值 也 能 够 一 定 程度 上 反映 出 GEO 的 近 圆 特 性 
经 大 量 仿真 和 实测 数据 验证 , 与 采用 固定 静止 
道 半 长 轴 标 称 值 相 比 , 该 优化 算法 的 定 轨 精度 更 
my, 且 更 加 稳健 . 


4 ”实验 

利用 位 于 郑州 的 工 型 正 交 双 基 线 CEI 测 轨 系 统 ， 
对 亚太 七 号 (Asia-Pacifc STAR 7, APSTAR 7) 进 
行 短 弧 定 轨 和 预报 实验 . CEI 的 东西 和 南北 基线 分 
别 为 75 m 和 35 m， 亚 太 七 号 下 行 Ku 信 标 信号 的 频 
点 为 12.442 GHz. 由 于 GEO 机 动 后 的 数 小 时 通常 不 
提供 精密 轨道 , 为 了 验证 算法 分 别 采用 了 CEI 的 仿 
真 数 据 和 实测 数据 . 


4.1 ”仿真 实验 

轨道 仿真 : 包含 东西 机 动 段 在 内 的 1 d 数 据 . 初 
始 时 刻 为 UTC (Coordinated Universal Time) 
2022-06-01 02:00:00, 机 动 时 段 为 当天 的 07:10:06 一 
07:30:06, 共 20 min. 机 动 采用 持续 小 推力 , 在 近 地 
点 附近 实施 . 仿真 采用 卫星 工具 包 (Satellite Tool 
Kit, STK) 的 Astrogator 模 块 , 顾及 了 10x10 阶 地 球 
引力 位 、 日 月 引力 摄 动 和 光 压 . 关键 时 刻 的 轨道 根 
数 见 表 1. 


法 分 


表 1 APSTAR 7 卫星 的 开 普 勒 轨道 根 数 
Table 1 Keplerian orbital elements of APSTAR 7 


Satellite Status Crucial Time (UTC) — a/km e i/? w/?^ Q/° Mf? 

Initial Orbit 2022-06-01 02:00:00 | 42165.9 0.00029 0.1506 357.2 82.0 276.5 

Start of East-west Maneuver 2022-06-01 07:10:06 | 42167.9 0.00036 0.1516 356.6 81.8 355.1 

End of East-west Maneuver 2022-06-01 07:30:06 | 42163.9 0.00027 0.1517 358.0 81.8 358.7 
图 3 是 轨道 半 长 轴 和 轨道 倾角 的 时 序 . 可 以 看 H, 采样 率 为 10 Hz. 假定 初始 模糊 度 已 知 , 且 暂 不 

出 , 非 机 动 期 间 , 半 长 轴 有 显著 的 半天 周期 震荡 , 机 考虑 模糊 度 失 锁 或 周 跳 现 象 . 

动 期 间 轨道 半 长 轴 快 速 下 降 了 4 km, 轨道 倾角 则 几 机 动 后 轨道 恢复 及 其 精度 评定 : 利用 机 动 结束 
平 无 变化 . 后 30 min 内 的 CEI 数 据 解 算 弧 段 中 间 时 刻 的 初 轨 改 
CEI 数 据 仿真 : 依据 CEI 的 相位 测量 方程 和 正 , 并 外 推 10 min. 仿真 轨道 作为 真 值 用 于 评估 短 
7.5 mm 随机 噪声 水 平 , 仿真 生成 1d 的 卫星 跟踪 数 弧 定 轨 - 预 报 段 的 轨 位 精度 . 需要 说 明 的 是 , 为 了 保 
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据 , 统计 轨道 解 算 精度 . 
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Fig.3 Semi-major axis and inclination series including the 


east-west maneuver 


在 无 轨道 先 验 约束 条 件 下 ，CEI 短 弧 定 轨 不 收 
SX. 依据 (10) 式 增加 了 部 分 参数 的 协 方差 约束 条 件 
后 短 弧 定 轨 均 收敛 . 

令 半 长 轴 的 初 值 在 静止 轨道 的 半 长 轴 标 称 值 
附近 一 20 km 一 20 km 窗口 内 变化 , 短 弧 定 轨 的 误差 
均 方 根 值 (Root Mean Square, RMS) 随 之 变动 ， 如 
图 4 所 示 . 图 中 R、T 和 N 分 别 是 初 轨 位 置 误差 的 径 
向 、 切 向 和 外 法 向 分 量 . 显然 , 短 弧 定 轨 的 精度 对 
半 长 轴 的 初 值 采 用 值 异常 敏感 , 尤其 是 R 位 置 分 量 ， 
这 与 CEI 的 测 角 性 质 一 致 . 此 外 ,可 以 看 到 初 值 取 
标 称 值 时 并 非 对 应 于 短 弧 定 轨 误 差 最 小 的 情况 . 

为 此 , 采用 自 适 应 优化 求解 半 长 轴 初 值 的 经 验 


黄 俊 迦 等 : 基于 CEI 测 量 的 GEO 


标 快速 轨道 恢复 


1 期 


证 蒙特 卡 洛 仿真 的 算法 稳定 性 , 仿真 了 20 次 CEI 数 


算法 . 图 5 给 出 了 初 值 在 半 长 轴 标 称 值 附 近 士 20 km 
范围 内 求解 的 偏心 率 时 序 . 令 偏 心率 取 极 小 值 时 
对 应 最 佳 的 半 长 轴 初 值 ， 可 见 该 值 比 标 称 值 低 约 
15 km. 
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Fig.4 RMS of the positional errors of the short-arc orbit 
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Fig.5 Initial values of semi-major axis and the solved 


eccentricities of the short-arc orbit determination 


分 别 采 用 


半 长 轴 的 标 称 值 和 优化 值 作 为 初 值 


进行 短 弧 定 轨 和 预报 . 表 2 为 20 次 仿真 解 算 的 轨 位 
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误差 RMS 统 计 ， 图 6 为 其 中 一 次 解 算 的 轨 位 误差 时 
序 , 图 中 实 线 和 虚线 分 别 对 应 30 minh Ag Jf x SUBE 
和 10 min 的 预报 段 . 可 以 看 


本 质 ，R 分 量 的 轨道 恢复 精度 始终 低 于 工 和 N 分 量 ， 
采用 固定 静止 轨道 半 长 轴 


标 称 值 


H: (1) 由 于 CEI 的 测 角 


时 ( 


图 a)，R 分 量 


误差 接近 10 km; (2) 采 用 自 适应 优化 的 半 长 轴 初 值 
E, R 分 量 改 善 最 为 显著 , R、T、N 位 置 分 量 误 差 


的 RMS 均 优 于 3 km; (3) 优 化 后 的 半 长 


动 性 不 大 ( 见 表 2). 


表 2 仿真 数据 的 短 弧 定 轨 和 预报 误差 RMS 统 计 


Table 2 Error RMS statistics of the short-arc orbit determination and orbit prediction with simulation 


的 计算 平均 值 比 标 称 值 小 14.76 km, 且 该 数值 的 波 


初 值 , 20 次 


Positional error / km 


data 
Initial Value Method [Initial Value of Semi- OD (Orbit Determination) and RMS/m 
of Semi-major Axis major Axis/km OP (Orbit Prediction) R T N 
OD 9174.49 | 739.93 1051.86 
Nominal Value 42165.76 
OP 9337.71 478.84 861.93 
OD 1784.15 183.03 714.18 
Optimized Value 42151.00 
OP 2069.15 157.46 1563.06 
ERE EN] 
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长 轴 初 值 对 应 的 短 弧 定 轨 和 预报 精度 


Fig.6 The short-arc orbit determination and orbit prediction accuracy with two adopted semi-major axes for simulation data 


4.2 ”实测 数据 实验 


目前 机 动力 的 高 精度 建 模 尚 有 困难 , 机 动 后 通 
常 需要 重新 累积 测 轨 数 据 , 因此 机 动 后 的 数 小 时 内 


TERM KIE SLPT?2. 为 了 模拟 机 动 后 的 快速 


轨道 恢复 场景 , 在 亚太 七 号 的 非 机 动 期 间 (2019 年 1 
月 7 日 至 10 日 ) 分 别 截取 了 4 段 CEI 数 据 进行 短 弧 定 
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轨 , 并 利用 国家 授时 中 心 提供 的 基于 测 距 的 长 弧 定 
轨 结 果 ( 米 级 精度 ) 进 行 外 符合 精度 评估 . CEI 测 量 


系统 的 配置 和 亚太 七 号 的 定点 经 度 等 与 仿真 实验 


一 致 . 


采用 相同 的 先 验 信息 约束 进行 30 min 短 弧 定 


轨 和 10 min 预 报 , 考察 了 两 种 半 长 轴 初 值 方案 : 
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方案 A 一 采用 静止 轨道 的 半 长 轴 标 称 值 ; 是 R 分 量 , 不 仅 是 短 弧 定 轨 段 , 还 极 大 改进 了 预报 精 


方案 B 一 采用 自 适应 优化 的 半 长 轴 . 
图 7 和 表 3 给 出 了 4 次 实验 的 短 弧 定 轨 - 预 报 精 
度 , 图 中 pos 表 示 位 置 误差 . 对 比 A、B 两 个 方案 , 可 
以 看 出 : (1) 自 适应 优化 得 到 的 半 长 轴 初 值 较为 稳 
定 ， 比 方案 A 的 固定 标 称 值 小 7.5-10 km; (2) 优 化 
后 的 半 长 轴 初 值 对 轨道 恢复 精度 改善 显著 , 尤其 
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图 7 实测 数据 的 短 弧 定 轨 和 预报 精度 比较 (A: K 


Fig.7 Comparisons of the short-arc orbit determination and 


静止 


RMS of the positional errors of the satellite / km 


RMS of the positional errors of the satellite / km 


度 , 三 维 位 置 整体 预报 精度 从 A 方案 的 5 km 减 小 到 
3.5 km 左右 , 前 3 段 甚 至 优 于 1.5 km. 

因此 , 对 于 非 合作 GEO 目 标 公 里 级 的 轨道 恢复 
精度 来 说 , 采用 自 适 应 优化 的 半 长 轴 初 值 的 B 方 案 
优 于 采用 静止 轨道 半 长 轴 标 称 值 的 A 方案 , 能 够 用 
于 机 动 后 的 快速 轨道 恢复 . 
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orbit prediction accuracy with real data (A: With nominal value, 


B: With optimized value) 
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表 3 实测 数据 的 短 缴 定 轨 和 预报 误差 RMS 统 计 


Table 3 Error RMS statistics of the short-arc orbit determination and orbit prediction with real data 


A Initial Value Method Initial Value of OD and RMS/m 
OS 
ii of Semi-major Axis  Semi-major Axis/km OP R T N p 
OD 4209.97 1659.86 2109.23 4992.78 
1 Nominal Value 42165.76 
OP 4928.94 245.87 1317.85 5108.00 
UTC 2019-01-07 
OD 766.24 1292.11 1649.12 2230.76 
11:00—11:30 Optimized Value 42155.96 
OP 110.64 299.19 854.21 911.83 
OD 2329.14 1136.21 1595.55 3043.30 
2 Nominal Value 42165.76 
OP 3139.50 130.47 433.28 3171.94 
UTC 2019-01-08 
OD 1541.17 859.99 1253.22 2164.56 
11:00—11:30 Optimized Value 42158.16 
OP 673.81 | 340.99 174.22 775.01 
OD 3412.46 1665.42 574.69 3840.41 
3 Nominal Value 42165.76 
OP 4241.84 535.08 631.52 4321.84 
UTC 2019-01-09 
OD 884.48 1336.88 411.74 1655.01 
11:00—11:30 Optimized Value 42157.36 
OP 128.72 216.83 1024.19 1054.79 
OD 2586.90 1057.04 1198.46 3040.67 
4 Nominal Value 42165.76 
OP 3415.97 2040.27 2045.97 4474.10 
UTC 2019-01-10 
OD 1484.14 1356.74 1566.29 2548.86 
11:00—11:30 Optimized Value 42157.76 
OP 599.35 2371.95 2426.58 3445.82 
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Rapid Orbit Recovery of the GEO Targets Based on CEI 
Measurement 


HUANG Jun-jia | DU Lan LIU Zejun ZHANG Zhong-kai ZHOU Pei-yuan LIU Long-di 


(College of Geospatial Information, Information Engineering University, Zhengzhou 450001) 


ABsTRACT Connected Element Interferometry (CEI) has been used for space surveillance as an all- 
day, all-weather passive tracking technique. Geostationary Earth Orbit (GEO) satellite requires frequent 
maneuver planning to maintain its subsatellite position or accomplish other tasks, which makes rapid orbit 
recovery extremely important for surveillance and early warning whenever a maneuver ends. This paper is 
to investigate CEI-based GEO short-arc orbit determination. Firstly, rank-deficiency in the conventional 
orbit determination algorithm is demonstrated and a priori orbit constraint is given to strengthen the 
positive definiteness of normal matrix. Then an adaptive optimization method to solve for the initial value 
of semi-major axis is proposed. Finally, two short-arc orbit determination and prediction experiments 
are carried out using both the simulated and measured CEI data of Asia-Pacific STAR 7 (APSTAR 7) 
satellite. The results show that an accuracy of 4 km is achievable for 30-min short-arc orbit determination 
and 10-min orbit prediction by using the proposed algorithm, which can meet the demands of rapid orbit 
recovery for non-cooperative GEO targets after maneuvers. 


Key words GEO (Geostationary Earth Orbit) satellites, non-cooperative objectives, CEI (Connected 
Element Interferometry), rapid orbit recovery 
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